December 12, 2002 
USC (NT) -02-2, SNUTP 02-014, TRI-PP-02-07 



Solar-neutrino reactions on deuteron 
in effective field theory 

S. Ando^t, Y. H. Song( 6 », T.-S. Park^, H. W. Fearing^, 

and K. KuboderaW H 

(a) Department of Physics and Astronomy, 
University of South Carolina, Columbia, SC 29208, U.S.A. 
^ Department of Physics, Seoul National University, Seoul 151-742, Korea 
< c > Theory Group, TRIUMF, Vancouver, BC V6T 2A3, Canada 



The cross sections for low-energy neutrino-deuteron reactions are calculated within 
heavy-baryon chiral perturbation theory employing cut-off regularization scheme. The 
transition operators are derived up to next-to-next-to-next-to-leading order in the Wein- 
berg counting rules, while the nuclear matrix elements are evaluated using the wave 
functions generated by a high-quality phenomenological NN potential. With the adop- 
tion of the axial-current-four-nucleon coupling constant fixed from the tritium beta decay 
data, our calculation is free from unknown low-energy constants. Our results exhibit a 
high degree of stability against different choices of the cutoff parameter, a feature which 
indicates that, apart from radiative corrections, the uncertainties in the calculated cross 
sections are less than 1 %. 



PACS : 25.30.Pt, 25.10. +s, 25.65. +t, 12.39.Fe 



t E-mail address : sando@nuc003.psc.sc.edu 
■^E-mail address : singer@phya.snu.ac.kr 
§E-mail address : tspark@nuc003.psc.sc.edu 
"E-mail address : fearing@triumf.ca 
II E-mail address : kubodera@sc.edu 



1 



1. Introduction 

This Letter is concerned with a theoretical estimation of the cross sections, a ud , for the 
neutrino-deuteron reactions 

v e + d — > e~ + p + p, v e + d^e + + n + n (CC), (1) 
vi + d^vi+p + n, ui + d^ui+p + n (NC), (2) 

where CC and NC stand for the charged-current and neutral-current reaction, respectively, 
and / denotes the lepton flavor (7 = e, /i, r). Recent SNO experiments [1, 2] have provided 
strong evidence for v e oscillations. In interpreting the existing and future SNO data, 
accurate estimates of a vd in the solar neutrino energy region {E v < 20 MeV) are of great 
importance. 

Recently, two theoretical approaches have been used for evaluating o~ vd . One is a tradi- 
tional method in which nuclear electroweak processes are described in terms of one-body 
impulse approximation (IA) operators and two-body exchange-current (EXC) operators 
acting on non-relativistic nuclear wave functions. The EXC contributions are derived 
from one-boson exchange diagrams [3], while the nuclear wave functions are obtained by 
solving the Schrodinger equation involving high-quality realistic nuclear interactions. For 
convenience, we refer to this method as the standard nuclear physics approach (SNPA). 
The successful applications of SNPA are well documented in the literature [4] . A detailed 
calculation of a vd based on SNPA was carried out by Nakamura, Sato, Gudkov and Ku- 
bodera (NSGK) [5], and this calculation has recently been updated by Nakamura et al. 
(NETAL) [6]. 1 

The second approach is based on effective field theory (EFT), which has been gaining 
ground as a new tool for describing low-energy phenomena in few-nucleon systems [9, 
10, 11]. Butler, Chen, and Kong (BCK) [12] applied EFT to the vd reactions, using the 
regularization scheme called the power divergence subtraction (PDS) [13]. Their results 
agree with those of NSGK in the following context. The EFT Lagrangian in PDS involves 
one unknown low-energy constant (LEC), denoted by L^, which represents the strength 
of axial-current-four-nucleon contact coupling. BCK adjusted L^a to optimize fit to the 
a ud of NSGK and found that, after this adjustment, the results of the EFT and SNPA 
calculations agree with each other within 1% over the entire solar- v energy region for all 
of the four reactions in Eqs.(l) and (2). Furthermore, the best-fit value of L ly ± was found 
to be of a reasonable magnitude consistent with the "naturalness" argument [12]. 

The fact that the results of an ab initio EFT calculation (with one free parameter 
fine-tuned) are consistent with those of SNPA is considered to give strong support for 
the basic soundness of SNPA. At the same time, it highlights the desirability of an EFT 
calculation of a vd free from an adjustable parameter. In this Letter we describe an attempt 
toward such a goal. We employ here a formalism recently developed in the studies of 
the solar Hep process and the solar pp fusion reaction [14, 15]. In this method, invoking 
heavy-baryon chiral perturbation theory (HB;\BT) , we construct transition operators from 
irreducible diagrams according to Weinberg's counting scheme [9]; the nuclear matrix 

1 For earlier calculations, sec, e.g., [7, 8]. 
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elements are evaluated by sandwiching the EFT-controlled transition operators between 
the nuclear wave functions that have been obtained by solving the Schrodinger equation 
involving high-quality realistic nuclear interactions. For convenience, we refer to this 
EFT- motivated approach as EFT*. It is known [14] that, for the present purposes, it is 
sufficient to consider up to next-to-next-to-next-to-leading order (N 3 LO) in HB^PT, and 
that to this order there is only one unknown LEC, denoted by cIr in [15]. Like L\a in [12], 
(Ir controls the strength of the axial-current-four-nucleon contact coupling and subsumes 
short-distance physics that has been integrated out. An important point noticed in [15] 
is that, since the tritium /5-decay rate rf is also sensitive to djt, we can determine d,R 
from the well-known experimental value of rf . Once d R is determined, we can make a 
parameter-free calculation of a udl and the purpose of this communication is to describe 
such a calculation. 2 We shall show that, apart from radiative corrections for which we 
refer to the literature [17, 18, 19], a vd given here is reliable with ~ 1 % precision. 

2. Calculational method 

For low-energy processes, we can work with the current-current interaction: 

H=^=Jd 3 x [V ud 4 cc \x)l^(x) + jj» c \x)l^(x)} , (3) 

where G' F = 1.1803 x 10 -5 [GeV~ 2 ][20] is the weak coupling constant, and V ud = 0.9746 is 
the K-M matrix element. G' F includes the inner radiative correction: G' F 2 = G 2 F (1 + A^), 
where Gf = 1.1166 x 10~ 5 [GeV~ 2 ] is the Fermi constant and is the inner radiative 
correction [20]. 3 The CC- and NC-lepton currents, /( cc ^ and /< NC ^, are well known; the 
CC- and NC hadronic currents, Jjf^ and Jj^ c \ are written as 

4 CC) (^) = V±{x)-A±{x), (4) 
J^ c \x) = (l-2sm 2 9 w )V°(x) - Al{x) -2sm 2 9 w V?(x), (5) 

where and represent the vector and axial current, respectively. The superscripts, ± 
and 0, are the isospin indices of the isovector current and S denotes the isoscalar current; 
9w is the Weinberg angle, sin 2 9w= 0.2312. 

The v-d reactions can lead to various values of the relative orbital angular momentum, 
L, of the final two nucleons. We concentrate here, however, on the L=0 state ( 1 5 , ), since 
it is this partial wave that involves the d R term and since the contributions of higher 
partial waves are well understood in terms of the one-body operators. The contributions 
from L > 1 are significant in the upper part of the solar neutrino energy region, 4 but their 
uncertainty is small enough to be ignored in the present context. 

The one-body (IB) currents can be obtained from the phenomenological form factors 

2 Similar parameter- free calculations have been carried out for the solar pp- fusion reaction and the 
solar hep process [15], and for p,-d capture [16]. 

3 For more detail discussion of the radiative correction, sec [6]. 

4 The L>1 contributions increase a vd by -3.8 % at E„= 20 MeV [6, 21]. 
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of the weak- nucleoli current. 5 The isovector vector and axial- vector currents are given in 
momentum space as 6 



9v{q)i^ + 9M(q) 



2,171 jv 



J a /(q) = u(p>) 



T 



9A(qh' I l5 + gp(q) 75 



u(p), 
u(p), 



(6) 
(7) 



is the isospin index, u(p) is the Dirac spinor for the nucleon, and m M (m^) 
is the muon (nucleon) mass; gv(q), <?m(<z), 9A(q), and gp(q) are the vector, magnetic, 
axial-vector, pseudoscalar form factors, respectively. It is known empirically that the 
first three form factors can be parametrized very well in the dipole form with the use 



of effective radii, r-y = 0.59 



M 



0.80 and r\ 



0.42 fm 2 , for g v (q), gui^q) and g A (q), 



respectively [27]. We are adopting here the usual normalization: gy(0) = 1, <?a(0)=<? j 4= 
1.267, and <7m(0)=«v = 3.706. Although gp(q) is not well known empirically, it is strongly 
constrained by chiral symmetry; an HB^PT calculation up to NNLO [22] leads to 



9p(q) = — 



mi 



-,g A m^m N r A , 



(8) 



where g n N = 13.5. In fact, the contribution of the gp term is tiny in our case. We apply 
a non-relativistic expansion of the above expressions and retain terms up to 0(l/m%) 
(corresponding to N 3 LO of the chiral order); the details will be described elsewhere [28]. 

The two-body (2B) current operators are derived from the chiral lagrangian £, which 
is expanded as £ = J2p £p — £o + C\ + • • where Co and C\ are LO and NLO lagrangians, 
respectively. Their explicit expressions are: 



C 



,V[/V-Dh 2/ (j \ S ■ A] A' J | i-A-A-^). 



(9) 



2m 



-N 



N 



'y . D f _ D i + 2 g A {v ■ A, S ■ D} - (8c 2 - g 2 A )(v ■ A) 2 



-8c 3 A • A - (4c 4 + 1)[^, S v ] [A M , A,] - 2i(l + k v )[S», S% 



N 



gA 



m N fl 



-MdiNS ■ ANNN + 2id 2 e abc e^ a pv^A a ' u NS a T b NNS (3 T c N} , (10) 



where v M is the velocity vector 



(1,0) and is the spin operator 2S 71 = (0, cr). 



The explicit expressions of the fields, D M , A^, /+,, and x+j are given in [16], and f w is 
the pion decay constant. The LEC's, q, have been determined by Bernard et al. at the 



5 The low-energy structure of the form factors has been studied in detail within HBxPT [22]. At N 2 LO, 
however, we in principle need to consider off-shell form factors [23], a feature that reflects arbitrariness in 
choosing fields [24]. The influence of off-shell terms, however, should be small at low energies; see, e.g., 
Ref. [25] . An EFT study of the one-body Gamow- Teller matrix element in the two-nucleon system [26] 
explicitly shows that the off-shell effects are sufficiently small for our present purposes. 

6 Since we consider the final 1 S'o state only, there is no contribution from the isoscalar current. 
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tree-level [29] ; 7 

c 2 = 1.67 ±0.09, c 3 = -3.66 ±0.08, c 4 = 2.11 ± 0.08. (11) 

The LEC's of the contact terms d 1:2 will be discussed later in the text. 

We construct 2B transition operators from 2B irreducible Feynman diagrams up to 
N 3 LO in Weinberg's counting rule [9]. Since the tree- level 2B operators are higher in 
chiral counting than the tree-level IB operators by two orders, we can limit ourselves to 
tree diagrams for the 2B operators. In addition, since the gp term is highly suppressed, 
we do not consider it in the 2B operators. 




Figure 1: Diagrams for two-body current operators of order v — 1 (a,b) and v=2 (c,d,e,f). 
The wavy lines with V and A attached denote the vector and axial-vector current, re- 
spectively, the dashed line denotes the pion, and vertices without (with) "X" arise from 
the LO (NLO) lagrangian. 

The diagrams for the 2B operators are given in Fig. 1. Since we only have nucleons 
and pions in C, the effects involving exchange of heavier mesons such as the a and p 
mesons are embedded in the contact term, diagram (f) in Fig. 1. We denote by A a 
momentum scale below which our nucleon-pion-only description is expected to be valid. 
To prevent the exchanged momentum from surpassing A, we introduce the cutoff function 
S\(k) = e~ k ^ 2A ^ in calculating the Fourier transforms of the 2B transition operators [15]. 
As noted in [15], the short-range part of the 2B contributions can be lumped together 
into an axial-current-four-nucleon contact coupling term with the strength d R , where 
d R = d\ + 2d 2 + |c 3 + |c 4 + |. Then, for a given value of A, we can determine d R from 
the empirical value of rf . The results are [15]: 

d R = 1.00 ±0.07, 1.78 ±0.08, 3.90 ±0.10, (12) 

for A = 500, 600, 800 MeV, respectively. The explicit expressions of the current operators 
for the CC reaction have been given in [16]. 8 



7 The relation between our dimcnsionlcss LEC's, fij's, and q's used in the literature is q = m^a. 
8 Insofar the final two-nucleon partial wave is limited to s-wave, one can use the same expression for 
the NC reaction (with an appropriate change in the coefficient of the vector current) 
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3. The total cross section 

The total cross sections are calculated using the nonrelativistic formula: 

2p 2 k' 2 



/E'+(k>-E„y)7(^ FiZ > E,) 3^ lTl2 > (13) 



spin 



with the energy conservation relation valid up to 1/mjy, 



m d + E v - E' - 2m 



N 



m N 



P 2 + -AK + k' 2 - 2E v k'y) 



0, 



(14) 



where E v (£") is the energy of the initial neutrino (final lepton), p is the magnitude of the 
relative three-momentum between the final two nucleons, k' is that of the outgoing lepton 
(k' = \k'\), and y is the cosine of the angle between the incoming and outgoing leptons 
(y = k u -k'). F(Z,E') is the Fermi function and is the deuteron mass. The transition 
matrix T is decomposed as T = T± B + T 2B: where T± B and T 2B are the contribution of the 
IB and 2B operators, respectively. These will be evaluated with the use of the Argonne 
V18 potential [30]. 

Since the calculation of T 1B is standard [28], we give here only the explicit expression 
for T 2B : 



]=T 2B = Pxlo^Xim d J dr ^F 6 u (r) 3l (qr/2)^^-u d (r) 

u 'o( r ) ( u d( r ) ~ V2w d (r)^ - u {r) (u' d (r) - V2w' d (rf)] jo{qr/2) 



+ F 7 



+ F 8 u (r)j (qr/2) 



yi\(r) 



2/ia(0 



w d (r) + F 9 u o (r)j (qr/2)y 0A (r)u d (r) 



— * — * 

+ F w u (r)j (qr/2)y 2A (r)w d (r) + F u u (r)5 A (r)u d (r) 



1/2 



+ F 12 / dyu (r)j (yqr) 



1/2 



w d {r) 
V2 



(15) 



where (5 = (G' F V ud f/2 for CC and (3 = G' F 2 /A for NC. t = 3 1 - <r 2 , with a { being the 
i-th nucleon spin operator; Xim d and Xo,o are the spin wave functions for the deuteron 
and the final two nucleons, respectively. The radial function uq corresponds to the final 
two-nucleon s-wave, while u d and w d are the s-wave and d-wave radial functions of the 
deuteron; j L (qr/2) is the spherical Bessel function; q^ is a momentum transfer between 
the currents, g M = k'^ — k^ and q — \q\. Furthermore, 



Fa = 



9a 
P 

J 7T 



2 v-Jq~ 



g A (l + Ky) 
2m N f 2 



-9a 



6m N f% 



J, F 8 



-9a 



V2m N f% 



iqx (qx J) - 
J, 



9A 



;Q- Jq, 



Fa = 



2g A m\ (c 3 + 2c 4 + \ 



J, 



(16) 



F\q — 



2y/2g A ml (ca - c 4 - j) ^ 



-* 2g A d R 

^11 — JKU , -T12 



6 



where J7 M is the lepton current in momentum space, and 



<*a('-)= l^r^Slik), yoA(r) = / ^e*>jS<*) (17) 



and y 1A = -r£y 0A (r), y 2 \ = l£yo\(r) ; Vo M (r) is obtained by exchanging the 



1 d 



pion mass m n to L = ym 2 + (1/4 — y 2 )^ in Eq. (17). In the above expression we have 
neglected the small terms proportional to q 2 . 

4. Numerical results and discussion 

As mentioned, we consider in this work the contribution from the final two-nucleon 
s-wave only. The corresponding total cross section is denoted by <J^°. Table 1 gives 
calculated in EFT* for the four reactions in Eqs.(l) and (2). 



E v 


vd — > e pp 


vd — > e + nn 


vd — * vnp 


vd — > vnp 


2 


0.004 











3 


0.047 





0.003 


0.003 


4 


0.158 





0.031 


0.031 


5 


0.348 


0.029 


0.096 


0.094 


6 


0.625 


0.120 


0.204 


0.198 


7 


0.996 


0.284 


0.357 


0.346 


8 


1.463 


0.525 


0.558 


0.538 


9 


2.030 


0.846 


0.808 


0.774 


10 


2.697 


1.247 


1.106 


1.054 


11 


3.468 


1.727 


1.455 


1.378 


12 


4.342 


2.286 


1.853 


1.746 


13 


5.321 


2.922 


2.302 


2.157 


14 


6.405 


3.633 


2.800 


2.610 


15 


7.596 


4.418 


3.349 


3.104 


16 


8.892 


5.274 


3.947 


3.638 


17 


10.29 


6.200 


4.594 


4.212 


18 


11.80 


7.194 


5.291 


4.824 


19 


13.41 


8.252 


6.036 


5.474 


20 


15.13 


9.374 


6.830 


6.161 



Table 1: The total cross section (in units of 10~ 42 cm 2 ) for the vd reaction leading 
to the final two-nucleon s-state. For each of the four reactions in Eqs.(l) and (2), 
calculated in EFT* is shown as a function of the incident neutrino energy E v [MeV] . For 
the cutoff parameter, A = 600 MeV has been used. 



The results in Table 1 correspond to the case with A = 600 MeV, and we now discuss 
the cutoff dependence. In Figs. 2 and 3 we plot the ratio, £ = <7ib+2b/°'ib, where <J\b+2B 
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represents obtained with both the IB and 2B currents included while <r 1B represents 
a ud° obtained with the IB current alone. Fig. 2 gives £ for CC [yd — > epp), while Fig. 3 
shows £ for NC [yd — > vnp). The three lines in each figure correspond to different choices 
of A. As can be seen from the figures, exhibits extremely small A dependence, with 
only 0.02 % changes over a wide range of physically reasonable values of A (A = 500 — 800 
MeV). 

We now briefly discuss estimation of higher chiral-order effects. The expansion pa- 
rameter here is Q/A, where Q is the pion mass or the typical external momentum 
scale Qexti and A is the chiral scale cutoff, A « 600 MeV. It is common to assume 
Qext ~ f^n, but Qext m our case is the incident neutrino energy E„, whose maximum 
value is E™ ax ~20 MeV; thus E™ ax /A ~ 0.03 < 0.23 ~ m n /A. The actual numerical 
behavior of the chiral expansion in the present case may be typified by the results for the 
CC reaction [yd — > e~pp) at E v = 20 MeV. As far as the IB operators are concerned, the 
contribution to a L=0 of the LO terms amounts to 88.5 %, while the corrections due to 
the NLO, N 2 LO and N 3 LO terms are 8.8 %, -0.5 % and ~ 0.001 %, respectively. As for 
the 2B operators, the N 2 LO terms give a ~ 0.3 % correction, whereas the N 3 LO terms 
give a ~ 2.9 % correction. Thus, the overall behavior is consistent with convergence with 
respect to the expansion parameter, m^/A; the rather conspicuous 2.9 % correction of the 
N 3 LO 2B terms is comparable to (m^/A) 3 ~ 1.2 %, while the other terms are decreasing 
faster (almost in powers of E™ ax / A). Therefore a possible measure of corrections due to 
N 4 LO or higher-order terms is 2.9 %x(m 7r /A) ~ 0.6 %. 

The convergence property, however, can in fact be better than this. Since in our ap- 
proach the overall strength, d R , of the 2B operator is adjusted to reproduce Fg, the bulk 
of higher order corrections have already been effectively taken into account. In particular, 
the chiral-symmetry breaking terms (proportional to m n ) give energy-independent con- 
tributions, which are essentially incorporated into the effective d R . The derivative terms 
acting on the wavefunctions or the two-body operators may pick up the pion mass scale, 
but their effects at the tritium /5-decay energy are again essentially subsumed in d R . The 
remaining pieces of higher-order contributions are independent effects, and hence they 
are likely to be controlled by the parameter E v j A rather than m^/A. From this viewpoint 
it seems reasonable to adopt 2.9 %x(E™ ax /A) ~ 0.1 % as a measure of the higher-order 
corrections. Another measure of convergence is obtained as follows. A tenet of a cutoff 
EFT (such as used here) demands that, provided an enough number of terms are included 
in chiral expansion, the calculational results should be independent of choices of the cutoff 
parameter A (within a reasonable range). Thus, the sensitivity of the calculated a L=0 to 
A serves as an indicator of the importance of the contributions of the neglected higher 
order terms. This sensitivity, however, has been found to be extremely small (0.02 % 
variation) in our case. 

Although the above discussion suggests that higher-order effects (N 4 LO or higher) 
are reassuringly small, we make a brief comment on three-body (3B) operators, which 
represent a particular class of higher-order contributions. It is known (see Table I of 
the last article in Ref. [15]) that, at N 4 LO, there is a contribution to the GT transition 
from the 3B-operator, which we denote here by O gt {3B). Obviously, although O gt (3B) 
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10 12 14 16 18 20 
E v [MeV] 



Figure 2: The ratio £ for CC defined in 
the text. The results for three different 
choices of A are plotted. The vertical 
bars represent changes in £ as d R is var- 
ied within a range allowed by the exist- 
ing experimental errors in rf ; the rep- 
resentative results obtained for A = 600 
MeV are shown for three values of E v . 




10 12 14 
E v [MeV] 



Figure 3: The ratio £ for NC defined in 
the text. See also the caption for Fig. 2. 



contributes to rf , it plays no role in the two-nucleon systems. At N 4 LO, therefore, in 
renormalizing d R with the use of rf , one would need to subtract the contribution of 
Ogt(3B). Formally speaking, our present treatment is free from this complication, since 
both the determination of d R and the calculation of a v ^ are carried out within N 3 LO. 
However, to the extent that d R adjusted to reproduce rf effectively includes higher order 
contributions, the above-mentioned subtraction is still needed. Although a full solution 
of this problem would require a systematic N 4 LO calculation, it is reasonable to expect 
that the contributions of the 3B operators, and hence the uncertainties due to them also, 
lie within the above-discussed overall range of higher-order effects. 

These considerations lead to the estimation that the corrections due to the N 4 LO or 
higher-order terms should be of the order of ~ 0.1 %. We also note that, within SNPA, the 
3B contribution to rf was calculated explicitly and found to be negligibly small compared 
with the leading 2B terms [31]. 

Figs. 2 and 3 also show the uncertainty in £ due to the finite precision with which d R 
can be fixed from rf . In fact, the largest uncertainty in our present calculation comes 
from this origin, and yet it only amounts to ~ 0.5 % ambiguity in £. Based on these 
observations, we consider it safe to conclude that cr^^'s calculated here are reliable at 
the ~ 1 % level. 

Comparison of our EFT* results with those of the latest SNPA calculation by NETAL 
[6] has already been described in [6]. We therefore only mention here that cr^° in table 
1 agrees with cr^° of NETAL within 1 % accuracy (see Table 4 in [6]). As discussed, 
to the chiral order we are concerned with, cr^p calculated in EFT* should agree with 
that obtained in SNPA. Therefore a u d (including all final partial waves) in EFT* can be 
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identified, within 1 % accuracy, with a ud given in NETAL[6]. 

There have been attempts to directly apply EFT to nuclear systems with mass number 
A> 3 [10, 32]. Here, "directly" means that the nuclear wave functions are obtained in 
the framework of EFT instead of using phenomenological potentials. It will be interesting 
to employ this "direct" EFT approach for determining dn (or L\a) from rf and use the 
resulting value of cIr for recalculating a u d- 

To summarize, we have carried out an EFT* calculation (up to N 3 LO) to estimate 
a^d , the cross sections of the vd reactions leading to the final two-nucleon s-wave state. 
Our results agree, within 1 % accuracy, with those of the most recent SNPA calculation 
reported in [6]. In addition, we have found that the calculated cr^ exhibits very small 
cut-off dependence (only ~ 0.02 % variation). The corrections due to higher chiral order 
terms are estimated to be of the order of ~ 0.1 %. The prime uncertainties in the 
calculated c^ stem from the experimental errors in rf ; this uncertainty, however, is less 
than ~ 0.5 %. We therefore conclude that, apart from the radiative corrections for which 
we refer to the literature, the uncertainties in the calculated (J^ are less than 1 %. 
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